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1. Introduction 

In this paper we consider the so-called crystallographic Radon transform 
(or crystallographic X-ray transform) and totally geodesic Radon transform 
on the group of rotations 50(3). As we show both of these transforms natu- 
rally appear in texture analysis, i.e. the analysis of preferred crystallographic 
orientation. Although we discuss only applications to texture analysis both 
transforms have other applications as well. 

In section 2 we start with motivations and applications. In sections 3 and 
4 we develop a general framework on compact Lie groups. In section 5 we 
give a detailed analysis of the totally geodesic Radon transform on SO (3). 
In section [6] we compare crystallographic Radon transform on 50(3) and 
Funk transform on 5 3 . In section [7] we show non-invertibility of the crys- 
tallographic transform. In section 8 we describe an exact reconstruction 
formula for bandlimited functions, which uses only a finite number of sam- 
ples of their Radon transform. Some auxiliary results for this section are 
collected in Appendix. 

The Radon transform on 50(3) has recently attracted attention of many 
mathematicians. In addition to articles, which will be mentioned in our 
paper later we also refer to 0, [E], [S], [HJ, |2U, 

2. Texture goniometry 

A first mathematical description of the inversion problem in texture anal- 
ysis was given in [B] and [7J. Let us recall the basics of texture analysis 
and texture goniometry (see [3] and [1]). Texture analysis is the analysis of 
the statistical distribution of orientations of crystals within a specimen of a 
polycrystalline material, which could be metals or rocks. A crystallographic 
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orientation is a set of crystal symmetrically equivalent rotations between an 
individual crystal and the specimen. 

The orientation distribution by volume requires a measure of the volume 
portion AV/V of a polycrystalline specimen of total volume V carrying 
crystal gains with orientations within a volume element AG C G of the 
subgroup G of all feasible G S 50(3). 

The orientation g of an individual crystal is the active rotation g £ 50(3) 
that maps a right-handed orthonormal co-ordinate system fixed to the speci- 
men onto another right-handed orthonormal co-ordinate system fixed to the 
crystal. In other words if a unit direction is represented by unit vectors h 
with respect to the crystal frame where as the same direction is represented 
by the vector r with respect to the specimen frame then h = gr. 

Assuming that the measure possesses a probability density function / : 
G ->■ M+, then 



prob( 5 G AG) = / f(g) dw 
Jag 



AG 

and / is referred to as the orientation density function (ODE) by volume 
and the Riemannian measure on 5 3 is 

dijjg = sin f3dad[3d'y, 

where a, (3, 7 are the Euler angles. Note that du g = 8ir 2 dg, where dg is the 
invariant Haar measure on 50(3). 

In X-ray diffraction experiments, the orientation density function / can- 
not be directly measured but with a texture goniometer the pole density 
function (PDF) P(h, r) can be sampled. P(h, r) represents that a fixed 
crystal direction h or its antipodal —h statistically coincides with the spec- 
imen direction r due to Friedel's law in crystallography [10]. The PDF or 
crystallographic X-ray transform is the tomographic projection of an orien- 
tation density function given by the formula 

(1) Pf(h, r) = \(Mf(h, r) + Mf{-h, r)), 
where the Radon transform is defined as 

(Mf)(h,r) = ±- f f(9)dw g = 

2vr J{geSO(3):h=gr} 

(2) 4n f f(g)5 r (g~ 1 h)dg = (f * 5 r )(h), (h, r) G 5 2 x 5 2 , 

JSO(3) 

where 5 r is the measure concentrated on the set of all g 6 50(3) such that 
h = gr. 

Since ODF / is a probability density it has to have the following proper- 
ties: 

(a) /( 5 )>0, h,reS 2 , 

( b ) Jso(3) f(9)dg = 1. 
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In what follows we will discuss inversion of crystallographic X-ray trans- 
form Pf and Radon transform Mf. 

First we formulate what can be called analytic reconstruction problems. 

Problem 1. Reconstruct the ODFf(g), g G SO(3), from PDF Pf(h,r), h,r G 
S 2 . 

It will be shown in section [7] that this problem is unsolvable in general 
since the mapping 

f^Pf 

has non-trivial kernel. 

Problem 2. Reconstruct f(g), g G SO(3), from all&f(h,r), h,r G S 2 . 

An explicit solution to this problem will be given in section 
In practice only a finite number of pole figures P(h,r), h,r G S 2 , can be 
measured. Therefore the real life reconstruction problem is the following. 

Problem 3. Using a finite number of pole figures P(hi,rj), h{,rj G S 2 , i = 
1, . . . , n, j = 1, . . . , m, find a function f on SO(3), which would satisfy (in 
some sense) equations and {]|) and conditions (a) and (b). 

An approximate solution to this problem, which is described in terms of 
Gabor frames was found in [8]. 

The corresponding discrete problem for Mf can be formulated as follows. 

Problem 4. Reconstruct f(g), g G £0(3), from a finite number of samples 
&f{hj,rj), hj,rj G S 2 , j = 1, ...,m. 

This problem will be solved in section [8] for bandlimited functions on 
SO (3). We were able to obtain an exact reconstruction formula for bandlim- 
ited functions, which uses only a finite number of samples of their Radon 
transform. 

In section H] we extend Problem 2 to general compact Lie groups. Namely, 
we consider the following problem. 

Problem 5. Let be a subgroup of the compact Lie group if. Determine 
domain and range for the Radon transform 

&f(x,y) = I /(xhy-^dh, x,ye<f. 

In section [3] we remind basic facts about Fourier analysis on compact 
Lie groups. In section [6] we compare crystallographic X-ray transform on 
SO(3) and Funk transform on S* 3 . In Appendix [9] we briefly explain major 
ingredients of the proof of our Discrete Inversion Formula in section [HJ 

3. Fourier Analysis on compact groups 

Let if be a compact Lie group. A unitary representation of if is a con- 
tinuous group homomorphism ir: if — > U(d 7T ) of if into the goup of unitary 
matrices of a certain dimension d n . Such a representation is irreducible if 
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n(g)M = Mir(g) for all g G ( S and some M G C d ^ xd ^ implies M = cl, 
where I is the identity matrix. Equivalently, C dw does not have non-trivial 
7r-invariant subspaces V C ( C dn with Tt{g)V C V for all g G . Two repre- 
sentations tti and 7T2 are equivalent, if there exists an invertible matrix M 
such that 7Ti(»M = Mn 2 (g) for all g G Sf. 

Let denote the set of all equivalence classes of irreducible representa- 
tions. This set parameterizes an orthogonal decomposition of the Hilbert 
space L 2 (&) constructed with respect to the normalized Haar measure. 

Theorem 1 (Peter- Weyl, [21] )• Let & be a compact Lie group. Then the 
following statements are true. 

a: Denote H n = {g H> trace(vr(g)M) : M G C^*^}. Then the Hilbert 
space L 2 (£f) decomposes into the orthogonal direct sum 

(3) L 2 (W) = ©F, 

b: For each irreducible representation ir G & the orthogonal projection 
L 2 (%?) — > H w is given by 

(4) f^d T f(h)xAh~ 1 g)dh = d n f*XiT, 

in terms of the character X-k{q) = trace("7r(g)) of the representation 
and dh is the normalized Haar measure. 

We will denote the matrix M in the equation / * Xir = trace(7r(g)M) as 
Fourier coefficient /(7r) of / at the irreducible representation tt. The Fourier 
coefficient can be calculated as 

/>)= / f(g)n*(g)dg. 
The inversion formula (the Fourier expansion) is then given by 
fid) = ^trace(7r(s0.f(7r)). 

If we denote by = trace(M*M) the Frobenius or Hilbert-Schmidt 

norm of a matrix M, then the following Parseval identity is true. 

Proposition 1 (Parseval identity). Let f G L 2 (^). Then the matrix-valued 
Fourier coefficients f G <C d7rXdjr satisfy 

(5) \\ftf=E d -\\fW\\HS- 

On the group one defines the convolution of two integrable functions 
/, r G L 1 ^) as 

f*r(g) = I f(h)r(h~ 1 g)dh. 
Since f*r G L 1 (^), the Fourier coefficients are well-defined and they satisfy 
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Proposition 2 (Convolution theorem on Let /, r G L 1 ^) then f *r G 
L^Sf) and 

/ * r(7r) = /(vr)r(7r). 

The group structure gives rise to left and right translations T s / i— >■ /(<? ■) 
and T 9 / i— > /(•#) of functions on the group. A simple computation shows 

fj(7r) = f(7r)7r*(g) and f^/(vr) = vr( 5 )/(vr). 

They are direct consequences of the definition of the Fourier transform. 

The Laplace-Beltrami operator Ag? of an invariant metric on the group 
<S is bi-invariant, i.e. commutes with all T g and T 9 . Therefore, all its 
eigenspaces are bi-invariant subspaces of L 2 (&). As H n are minimal bi- 
invariant subspaces, each of them has to be eigenspace of Ay with cor- 
responding eigenvalue — A^. Hence, we obtain 

A^/ = - ^2 A^trace(vr(5)/(7r)). 

4. Problem 5: Radon transform on compact groups 
4.1. Radon transform. 

Definition 1. Let be a subgroup of the compact Lie group <S . The Radon 
transform 3% of an integrable function f on is defined by 

(6) £f{x,y)= [ /(xhy-^dh, x,yeW, 

where dh denotes the normalized Haar measure on J^. 

The range of fflf should be contained in & x . 

Lemma 2 ([2]). The Radon transform (fjty is invariant under right shifts of 
x and y, hence the range is a subset of^/Jf? x tf/Jf?. 

Proof. First, we take the Fourier transform oiMf with respect to the x and 
let y be fixed and regard 3&f(x, y) as a function of x G 5f only then 

WTy)(T) = TT^TT* («)/>), 7T G <$. 

It is easily seen that &f(x, y) is invariant under the projection and we 
obtain 

&f(x-h,y)=£f(x,y) VheJlf. 

If we look at the Radon transform as a function in y while the first argument 
x is fixed, we find 

P*K^/)(x,y)= / <%f(x,yh)dh= f Vd 7r trace(/( 7 r)7r( ;C )) 7 r jr 7r(/ l - 1 y- 1 

tt£G 

= d 7r trace(/(7r)7r(a;))7r^7r*(y) = {M)f(x, y). 

7T£G 
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Consequently, 3&f{x, y) is constant over fibers of the form yj^ and 

□ 

We have just proven that this Radon transform act as a projection from 
5f — >• &/J$?. To fix the mapping properties we need function spaces. We 
start with L 2 (W/J^ x and get the following theorem. 

Theorem 3 ([2]). Let J4f be a subgroup of ' which determines the Radon 
transform on & and let <S\ C G be the set of irreducible representations with 
respect to Jff. Then for f € C°°(£f) we have 

(ttjt) \\f\\ 2 HS- 

Proof. We expand &f(x,y) for fixed y into a series with respect to x and 
apply Parseval's theorem 

I Wlliaw/jrx^/jf) = X^m \\ 7r ^ 7r *(y)f( 7r )\\ 2 Hs d y 

= dn J^ trace (f*( 7r ) 7T (y) 7T ^ 7r *(y)f( 7T ))dy 

= d 7r trace(/*(7r) ^ ^{y)^^* (y) dyj /(vr)) 

7T6Sf 

/ rank vr^ _^ \ d?r 

= ^ d 7r trace(/*(7r) I ^ J^n ik {y)^l. j y)dy J /( 



x 7T 

i,i=l 



2 ^trace(/*(7r)^^Id/>)) 

7r€3f 

^ rankTT^trace (f*(ir)f(ir)) 
J2 rank (tt^) ||/||^. 



□ 



From that we can conclude that the Radon transform is an isometry 
between L(£f) and some space on j .J? x 1 3^ but that doesn't say what 
is the right space. For that we start from the other side, i.e. we start with 
the L 2 -norm of /. 
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4.2. Mapping properties. We will need some help from harmonic analy- 
sis. The goal here is to find a Hilbert space structure such that the Radon 
transform is an isometry, which means that 

\f\\h(<$) = J2 d T I I/O) I Iks = d *\\f(' !r )\\h(&/jg'x&/j%') 

= lll^/lllla^/jrx^/jr) 
Suppose that the compact Lie group ^ acts transitively on a compact Rie- 
mannian manifold Ai as a group of isometries. Then if we set 

l(g)f(x) = fig-'x), ge&, x e M, 

we have a unitary representation of on L 2 (M). Now if A is the Laplace 
operator, we have 

l(g)A = Al(g), g£&. 

Thus the eigenspaces 

V x = {u € C°°{M) : Au = -X 2 u} 
of the Laplace operator, we have 

Kg) ■■ v x ->• v x , 

for each eigenvalue —A 2 . 

Lemma 4 (|25], page 120). If M is a compact rank one symmetric space, 
then acts irreducibly on each eigenspace Vx of A on M. 

Examples of such compact rank one symmetric spaces following |26j are 

& = SO(n+l), M = S n , 

^ = SU(n + 1), M. = CP n (complex projective plane ). 

Definition 2. A representation 1(g) is called a representation of class- 1 
relative to Jf? if in its space there are nonzero vectors invariant relative to 
J4? and the restriction of 1(g) to r Jt? is unitary. 

Lemma 5 ([26]. page 121, 122). If A4 = @/Jt? is a rank one symmetric 
space, with J^f connected, than I? (M) contains each class- 1 representation, 
exactly once, as an eigenspace of A. 

Therefore we will try to build up a norm from the eigenspaces of A which 
span some subspace of L 2 (M) and consider the case = SO(n + 1), Jf? = 
SO(n) in more detail. 

4.3. The case = SO(n + 1) and = SO(n). The following result is 
well-known. 

Lemma 6. The space 

L 2 (S n - 1 ) = ®H k , 

k 

where Hk is the space of harmonic polynomials, homogeneous of degree 
k, k£ N . 
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Both previously mentioned lemmas together are only fulfilled for S 2 , in 
this case we have 

Lemma 7 ([26j. page 140). The decomposition 

L\S 2 ) = §V k 

k 

contains each irreducible unitary representation of SO(3), exactly once. 

But also for higher dimensional spheres we have some nice results. We 
start with the orthonormal system of spherical harmonics yi G C°°(S n ), k G 
No, i = 1, • • • , dk(n) normalized with respect to the Lebesgue measure on S n . 
Obviously TLk = span{3^}f=j[ • Then the Wigner polynomials on SO(n + \) 
Tk {9)1 9 £ SO{n + 1) are given by 

r k \g) = I yi{g- 1 x)y{(x)dx 

JS" 

and due to orthogonality of the spherical harmonics 

d k (n) 

yl(g- l x)= £ Ti\g)y{(x). 
i=i 

From these properties and the orthonormality of the spherical harmonics 
it easy to see that the Wigner polynomials build an orthonormal system in 
L 2 (SO(n+l)). Unfortunately, Wigner polynomials do not give all irreducible 
unitary representations of SO(n + 1) if n > 2. (Compare with Lemma 0) 

Definition 3. Let £ be a linear space of functions f(x), f x G M. If 
in the space C of any representation of class- 1 relative there is only one 
normalized invariant vector, then J$? is called a massive subgroup. 

Lemma 8 ([27], Chapter IX. 2). SO(n) is a massive subgroup of SO(n + 
1). Furthermore, the family Tk, k G No, gives up to equivalence all class- 1 
representations of SO{n + 1) with respect to SO{n). 

For the following let xq be the base point of SO(n+l) / SO(n) ~ S n (xo is 
usually chosen to be the "north pole".) In this case the set of zonal spherical 
harmonics is one-dimensional and spanned by the Gegenbauer polynomials 

c (n-l)/2^T x )_ We recall 

some helpful and well known results. 
Lemma 9 (Addition theorem). For allx, y G S n , k G No andi = 1, . . . , dk{n) 

r (n-l)/2 f T x , Qn , d k {n) 

Lemma 10 (Zonal averaging). 

JsO(n) Cf >' (1) 
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Lemma 11 (Funk-Hecke formula). Let f : [— 1, 1] — )■ C be continuous. Then 
for all i = 1, . . . , dk{n) 

f(x T y)yi( X ) dx = yi(y) c l S m {1) £ mc { r iy \t)(i - w- 1 &. 

We look for all irreducible representations of SO(n + 1) which do not 
have vanishing coefficients under the projection Pso(n)- These are the class- 
1 representations of SO(n + 1) with respect to SO(n) and the projections 
are given by 

^Tt j =f V k j {g)dg=f f yi(g^x)dgy{(x)dx 

JSO(n) JS n JSO(n) 

f ct- l), \xlx)yi{x)dx 



c (n-l)/2 J s „ 



'k 

(cp 7 ^)) 2 
\s n 



yi(x )yi(x ) 



d k {n)' 

due to the Funk-Hecke formula and the normalization of Gegenbauer poly- 
nomials. We asume that the basis of spherical harmonics y k (x)is chosen 

such that y k (xo) = J d , k S™) and y k (xo) = for all i > 0, then 



, rn y *W = O^ 7 * )(*) = / r k \ g h) dh = T k a ( g ), x = gXo . 

V d k [n) JSO(n) 

Now we calculate the Radon transform from the span (T k ) k £^ into S n x S n 
Let be 

oo d k (n) 

M = E E 

k=0 i,j=l 

then 

oo 

@f(x,y) = ^2d k (n)trace(f(k)Tk(x) T k *(y)) 

k=0 

oo d k (n) 

k=0 i,j=l 
oo I ™j | d k (n) 



Et^4(«) E f(khyi(x)yi(y) 

k=Q a kVL) . j=1 



|S n lE E /Wy^Wftv) 

fc=0 i,j'=l 
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From the last equation it is easily seen that 

(7) A X M = A y 3!. 

Which also means that the image is a subspace of L 2 (S n x 5 n ) rather than 
the whole space. It also follows that for the Fourier coefficients of g = 3%f 
we have 

(8) g(k)a = \s n \f(k) l3 <=> f(k)ij = rLgWij- 

5. Problem 2: Radon transform on 50(3) 

Problem 3 can be derived from Problem 5 by choosing £f = 50 (3), ffl = 
50(2) and thus Sf/JT = 50(3)/50(2) = 5 2 . In this special situation all 
irreducible representations are equivalent to an irreducible component of the 
left regular representation 

T(g) : f(x)^f(g- 1 -x), 

where • denotes the canonical action of 50(3) on 5 2 . The T invariant sub- 
spaces of L 2 (5 2 ) are %k = {3^, i = 1, • • • , 2k + 1}, which are spanned by all 
spherical harmonics of degree k. It is well known that the eigenvalues of the 
Laplacian on 50(3) and on 5 2 corresponding to polynomials of degree k is 
-k(k + 1), i.e. 

A so{3) 7t j = -k(k + l)7t j and A s2 yi = -k(k + l)yl 

We will use that to define an appropriate norm for the Radon transform from 
L 2 (50(3)) onto an subspace of L 2 (5 2 x 5 2 ). Starting with the dimension 
of the representations space d n = 2k + I and —A? = — k(k + 1) we get 

4 = «/l + 4Ag and ^4 = ^2(A| + \\ + 1). Using © and the fact that A 
is equal to — k(k + 1) on the eigenspace Hu we obtain 

oo 

ll/lli 2 (50(3))=E(^ + 1 )H/>)llH5 = 
fc=l 

oo 

^(2m)||(4^)- 1 />)ll| 2 (5 2 x 52 ) = ll(^)" 1 (-2(A 1 +A 2 )+l) 1 /^ / ||| 2(52x52) , 

k=l 

where Ai + A 2 is a Laplace operator on 5 2 x 5 2 . Thus we define the following 
norm for u £ 0°°(5 2 x 5 2 ) 

IIMII 2 = ((-2A 5 2 x5 2 + l) 1/2 u, n) L 2 (5 2 x5 2), 

where A52 X 52 = Ai + A 2 . 

Because we have an isometry (up to a constant) the inverse operator is 
given by the adjoint of £%. To calculate the adjoint operator we express the 
Radon transform in another way. Going back to our problem in crystallog- 
raphy we first state that the great circle C XtV is given by 

C x>v = x'S0(2)(y 1 )- 1 := {x'h{y')-\ h G 50(2)}, x' , y' G 50(3), 
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and x' -xq = x, y' -xq = y with SO{2) being the stabilizer of xq G S 2 . Hence, 
&f(x,y)= f f(xhy- 1 )dh = 47r [ f(g)dg 

JSO(2) JC X ,y 

= 4tt / f{g)5 y {g- 1 • x) dg, f G L 2 (SO(3)), 

JSO(3) 

where C x>y = {g G SO (3) : g ■ x = y} is a great circle in 5*0(3) . To calculate 
the adjoint operator we will use the last representation of We have 

{M*U, f) L 2 {SO {3)) = ((-2A 5 2 XS 2 + l) l ' 2 U, <%f) L 2 {S 2 xS 2 } = 

(4tt) / {-2A S 2 XS 2 + I)u(x,y) [ f '(g)5 y (g- 1 ■ x) dg dx 
JS 2 xS 2 JSO(3) 

(An) f / (-2A S 2 xS 2 + l) 1 / 2 u(g-y, y)dyf(g)dg, 
JSO(3) JS 2 

i.e. the L 2 -adjoint operator is given by 

(9) M*u = (4tt) / (-2A s2xS2 + if'Mg ■ y, y) dy. 

Js 2 

Definition 4 (Sobolev spaces on S 2 x S 2 ). The Sobolev space Ht(S 2 x 

t 

S ), i 6 I, is defined as the domain of the operator (1 — 2A 5 2 X 52)2 with 
graph norm 

||/|| t = ||(l-2A 52xS2 )I/|| i2(s2xS2) , 

and the Sobolev space H^(S 2 x S 2 ), t G M, z's defined as the subspace of all 
functions f G #t(S 2 x S 2 ) suc/i A 1 f = A 2 f. 

and 

Definition 5 (Sobolev spaces on SO (3)). The Sobolev space Ht(S0(3)), t G 

t 

M, is defined as the domain of the operator (1 — 4Aso( 3 ))2 with graph norm 
\\\f\\\ t = ||(1 - 4A g0( 3 ) )I/|| i2(50(3)) , / g L 2 (SO(3)). 

Then: 

Theorem 12. For any t > the Radon transform on SO(3) is an invertible 
mapping 

(10) @ : H t (SO(3)) -> fl^i (S 2 x S 2 ). 

t+ 2 

and 

(11) /(<?)=/ (-2A + l)Haf)(gy,y)dy = ]-(&*af)( S ). 

JS 2 47T 

Proof. For the mapping properties it is sufficient to consider case t = 0. 
Because the Radon transform is only an isometry up to 4tt, we obtain (jlip . 

□ 
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Since 

<%(T k )(x,y)=T k (^ SO (2) (T k 
we have 

Aqr : 

#7*0*, v) = 73(x)7j(y) = ^rKW^fe). 

Theorem 13 (Reconstruction formula). Let 

oo 2k+l 

G(x,y)=af(x,y) = E E G{k) i3 yi{x)y{{y) G H$ +t {S 2 x 5 2 ), t > 0, 
fe=0ij'=l 2 

6e i/ie result of a Radon transform. Then the pre-image f G Ht(SO(3)), t > 
0, is giwen fry 

oo 2k+l /rti I -i \ 00 2fc+l 

/ = E E = E E( 2fc + i)/(*)«t3 

fc=0 i,J=l fc=0 i,j=l 

oo 

= J2( 2k + l)trace(/(A;)r fc ). 

fc=0 

6. The group 50(3) and the sphere 5 3 

bmce 5 3 is a double cover of 50(3) every function / on 50(3) can be 
identified with an even function on 5 3 . After such identification the crystal- 
lographic Radon transform on 50(3) becomes the Funk transform on 5 3 . 
In this section we show how some known results for Funk transform can be 
interpreted in crystallographic terms. 

Let C denote the set of all 1— dimensional totally geodesic submanifolds 
C C 5 3 . Each C G C is a 1-sphere, i.e. a circle with centre O. Each circle 
is characterized by a unique pair of unit vectors (h,r) G 5 2 x 5 2 . 

Thus, the pole density function can be rewritten as 



Pf(h,r) = ±- { f(q)dut(q), 
07T Jcuc 1 - 



where the circle C C 5 3 represents all rotations mapping h G 5 2 on r G 
5 2 and where C 1 - ((ft, 52) is the orthogonal circle representing all rotations 
mapping —h on r, and where wi denotes the usual one-dimensional circular 
Riemann measure. 

Following Helgason (see [12], [13]) [H])) we introduce the next definition. 

Definition 6. For a continuous function f defined on 5 3 its 1-dimensional 
spherical (totally geodesic) Radon transform IZf is a function, which is de- 
fined on any 1-dimensional totally geodesic submanifold C C 5 3 by the fol- 
lowing formula 

W){C) = ^J c f(q) <M<?) = j f(q) dm(q), 

with the normalized measure m = tt-cji . 
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The Radon transform of / may be represented as the convolution of / 
with the indicator function of the great circle C. It associates with the 
function / its mean values over great circles C € C. Since each great circle is 
uniquely characterized by a pair (h, r) £ S 2 xS 2 (see [21]) and its antipodally 
symmetric, we also use the notation (lZf)(h,r) = (1Zf)(—h,—r) whenever 
it is more instructive. 

It should be noted that no distinction has been made whether / refers 
to 50(3) or S s , even though the form of / depends on the representation 
of g; in particular, with respect to 5 3 only even functions / could be ori- 
entation probability density functions as q G S 3 and —q represent the same 
orientation. 

Thus, one has 

(12) Pf(h,r) = \((Kf)(C h , r ) + (tt/)(C_fc, r )), 

where P(h, r) is again the crystallographic X-ray transform. 

The following definition of the generalized 1-dimensional spherical Radon 
transform and the respective dual can be found in [13], [14j . 

Definition 7. The generalized 1-dimensional spherical Radon transform of 
a real function f : S 3 h-> Mi is defined as 

(^ p) f){C) = -JL— f f{q)dq. 

47T 4 Sinp Jd{q,C)=p 

The next definition was introduced in PQ, and [9], p. 64. 

Definition 8. The spherically generalized translation of a real function F 6 
C{S 2 ) or F G LP(S 2 ), 1 < p < oo, is defined as 

(Ti p) F)(r) = 1 [ F(r')dr'. 

Z7TV 1 — COS z p Jr-r'=cosp 

where 2-Ky/l — cos 2 p is the length of the circle c(r; p) centered at r with 
radius cos p. 

To describe connection between Radon transform and other operators it 
is useful to introduce the angle density function 

((Af)(h,r;p)=R^f)(C h , r ). 

Note, that (Af)(h,r; p) has been introduced in [6], p. 44 and [7], p. 74, but 
with a false normalization. 

The following statement holds. 

Proposition 3. The generalized 1-dimensional spherical Radon transform 
is equal to the translated spherical Radon transform and it can be identified 
with the angle density function 

(13) (T 2 ip) [Kf])(h,r) = (TZ^f)(C h , r ) = (Af)(h,r;p) 
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The proof can be found in [3] . 



The following properties should be noted 

(14) (Af)(h,r;0) = (Kf)(h,r), 

(Af)(h,r;ir) = (Kf)(h,-r). 

According to its definition the quantity (Af)(h, r; ir) it is the mean value of 
the spherical pole probability density function over any small circle centered 
at r. Thus, it is the probability density that the crystallographic direction 
h statistically encloses the angle p, < p < it, with the specimen direction r 
given the orientation probability density function /. Its central role for the 
inverse Radon transform was recognized in [T5] [19] . 

The next two definitions were given by Helgason in [T3], [14] . 

Definition 9. The dual 1 -dimensional spherical Radon transform of a real 
continuous function p : C \— > Mi is defined as 



{n^){ q )= / <p(C)d»(C), 

JCSq 

where p denotes the unique measure on the compact space C £ C : q £ C, 
invariant under all rotations around q, and having total measure 1. 

Thus (TZip) (q) is the mean value of ip over the set of circles C G C passing 
through q. 

Definition 10. The dual generalized 1-dimensional spherical Radon trans- 
form of a real continuous function (p : C i— > Ri is defined as 

(7^V)(<Z)= / p{C)dp{C). 

With the normalized measure Ji = A~ 1 (p)dp, A{p) = 47rsin 2 /3, it is the 
mean value of p over the set of all circles C £ C with distance p from q. 

One has that (lZ^ (TZf)j(q) is the mean value of (TZf) over the set of 

all circles C G C with distance p from q, i.e. tangential to the sphere 
s(q; p). More specifically, with the usual two-dimensional spherical Riemann 
measure L02, 

(£W (TZf)) (q) = ( (Kf)(C h>r ) dftC h>r ) 

V ' Jd(q,C)=p 



1 



47TSm(2/)) J S 2 Jc(gh;2p) 

(15) =\f (Af){h,gh-2p)duj 2 (h) 

= \ I (KMf)(C h , gh )dw 2 (h) 
2 Js 2 

(16) = (K(TZMf))(q), 



(lZf)(C Kr )drduo 2 {h) 
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which is instrumental for the inversion of the spherical Radon transforms. 

Now, we use the inversion formula of S. Helgason [T3], which was obtained 
for the general case two-point homogeneous spaces. For two dimensional 
sphere the totally geodesic Radon transform is also known as Funk trans- 
form. We start with the inversion formula and derive a connection to the 
angle density function Af. The inversion formula reads as 

Moreover, one can use angle density function Af to deduce Matthies' inver- 
sion formula, which was obtained in |19j). One starts with (j 151) to obtain 



1 



d f u f 

— / / (Af)(h,gh;2aTccosv)dw 2 (h)—^. 
du Jo Js 2 V« 2 - v 2 



dv 



(Af)(h,gh-Tv)du 2 (h) 

47T \J S 2 



+ 

Using the formula 



/ / — (Af)(h, gh; 2 arccos \/s — w) dh—= dw 
Jo Js 2 ds yjw 



j" (AO = W) 

as aw 



and taking into account 

(Af)(h,gh-7r) = (nf)(h,gh) 

(cf. Eq. p4p ). one obtains for 2w = 1 — cos 6 and s = 1 the desired inversion 
formula: 

f(q) = ^(J s2 (nf)(h,gh)du; 2 (h) 

Jo Js 2 a cost/ z / 

The practical importance of this formula is that Af is easily experimen- 
tally accessible and might yield an improved inversion compared to inverting 
just pole intensities. This formula was already obtained in [3] and is equiv- 
alent to the formula in [2D] which was given without any proof. 

7. Problem 1: On inversion of crystallographic X-ray 

transform 

Unfortunately, neither with the Radon transform M over SO (3) nor with 
the Radon transform 1Z over S s we are solving the crystallographic problem. 
The point is that since 

y k (- x ) = (-i) h yi(x), 
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one has for G(x, y) = &f{x, y): 

Pf(x, y)= l - (<%f(x, y) + &f(-x, y)) = \{G{x, y) + G(-x, y)) 



1 



oo 2fc+l oo 2fc+l 



2 E E G(fc)y34(^(i/) + E E G{k) i3 y k {-x)y{{ y ) 

\k=0i,j=l k=0i,j=l 
-. / oo 2fc+l oo 2fc+l 

= o E E G(k) t3 yi(x)yi(y) + E ^(%3^(^(y) 

\k=0i,j=l k=0 i,j=l 

oo 4Z+1 

= E E GWkyii(*)yii(y)- 

1=0 i,j=l 

In other words we loose half of the data needed for the reconstruction, 
because the experiment (which is measuring PDF Pf) only gives the even 
coefficients G(2l)ij. 

Since odd Fourier coefficients G(2l+l)ij of the function G(x, y) = &f(x, y) 
disappear one cannot reconstruct the function f(g),g £ SO(3), from Pf(x, y). 
Note that we have two additional conditions stemming from the fact that / 
is a probability distribution function: 



(1) f(g) > 0, i.e. / is non-negative, 

(2) Iso(3)f(9)dg = l. 

The second condition is just a normalization, the first condition is less 
trivial. 

We obviously can reconstruct the even part f e (g) from the even coefficients 
G(2l)ij. A legitimate question to ask is the following: Does there exist an 
odd part f (g) such that f(g) = f e (g) + f a (g) > and f a (g) is uniquely 
determined? 

8. Problem 4: Exact reconstruction of bandlimited functions 
on so (3) from a finite number of samples of their radon 

transform 

It is clear that in practice one has to face situations described in Problems 
3 and 4. Concerning Problem 3 we refer to [8] where approximate inverse 
was found using language of Gabor frames. A solution to the Problem 4 will 
be described in the present section. 

Let B((x,y),r) be a metric ball on S 2 x S 2 whose center is (x,y) and 
radius is r. 

As it explained in Appendix there exists a natural number Ng2 xS 2, such 
that for any sufficiently small p > there exists a set of points {(x u , y u )} C 
S 2 x S 2 such that: 

(1) the balls B((x u , y u ), p/4) are disjoint, 

(2) the balls B{{x v , y u ), p/2) form a cover of S 2 x S 2 , 



RADON TRANSFORMS ON SO (3) 17 

(3) the multiplicity of the cover by balls B((x u ,y u ), p) is not greater 
than N S 2 xS 2. 

Any set of points, which has properties (l)-(3) will be called a metric 
p- lattice. 

For an oj > let us consider the space E w (<SO(3)) of oj-bandlimited func- 
tions on SO(3) i.e. the span of all Wigner functions with k(k + 1) < u. 

The goal of this section is to prove the following discrete reconstruction 
formula (|18|) for functions / in E u (SO(3)), which uses only a finite number 
of samples of Mf . 

Theorem 14. (Discrete Inversion Formula) 

There exists a C > such that for any oj > 0, if 

then for any p-lattice {(x u , y u )}™=\ of S 2 x S 2 , there exist positive weights 

fJL„ X U~ 2 , 

such that for every function f mE w (50(3)) the Fourier coefficients c£ . (&f) 
of its Radon transform, i.e. 

Mf{ x , y) = J2 c h w) yi(x)yl(y)> K*> + 1) < ^ (*> v) e 52 x ^ 

are given by the formulas 

(17) clj [Mf) =J2p u {Mf) (x„, y v )yi(x u )y{{y u ). 

Moreover, the function itself can be reconstructed by means of the formula 
2k+1 (Ik 4-11 

(18) f(g) = E E J ^r lc t3 WWfo). 9 € 50(3), 

k i,j 

in which k runs over all natural numbers such that k(k + 1) < uj. 
Proof. As the formulas 

(19) Aso^rjf = -k(k + i)T^, A s2 yi = -k(k + i)yi 

and 



(20) ^(x,y) = w - r[ yi(x)yi(y) 

show the Radon transform of a function / G E CJ (5'0(3)) is w-bandlimited 
on S 2 x S 2 in the sense that its Fourier expansion involves only functions 

y%yl which are eigenfunctions of A S 2 x5 2 with eigenvalue — k(k + 1). Let 
E U {S 2 x S 2 ) be the span of Thus 

Tl : E„(SO(3)) -> «? W (S 2 x S 2 ). 
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Let {(xi,yi), (x m , y m )} be a set of pairs of points in SO (3) and M. v = 
x u SO(2)y~ 1 are corresponding submanifolds of 50(3), v = 1, ...,m. 

For a function / S E w (50(3)) and a vector (of measurements) v = (v^)™ 
where 




one has 

We are going to find exact formulas for all Fourier coefficients of fflf G 
£ui(S 2 x S 2 ) in terms of a finite set of measurement. Since SO (3) has di- 
mension three then according to Theorem [TU] (see Appendix) every product 

(Mf) y i k yl, where k(k + l) < u belongs to £ n (S 2 x S 2 ), where n = 4 x 6a; = 
24a; . 

By the Theorem [TJ] (see Appendix) there exists a positive constant C, such 
that if p = C(oj + l) -1 / 2 , then for any p-lattice {(xi, yi), (a; mw > J/m w > )} i n 
5 2 x 5 2 there exist a set of positive weights 

\i v x VT 2 

such that 

<fe W) = / W) y)yi(x)yi(y)dxdy = 

JS 2 xS 2 

N 

(21) E ^ w) y^M^M- 
i/=i 

Thus, 

V 

where 

TV 

(22) W) = 2 ^ W) (x u )yi(x u )yi( Xu ). 

v=l 

Now the reconstruction formula of Theorem [T3] gives our result (|18|) . 

□ 



9. Appendix 

We explain Theorems [TS] and [To! which played key role in the proof of 
Theorem Q3J 
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9.1. Positive cubature formulas on compact manifolds. We consider 
a compact connected Riemannian manifold M. Let B(£,r) be a metric ball 
on M whose center is £ and radius is r. 

It was shown in [23J, [24J, that if M is compact then there exists a natural 
number N^i, such that for any sufficiently small p > there exists a set of 
points {£fc} such that: 

(1) the balls J3(£fc,p/4) are disjoint, 

(2) the balls B(l k ,p/2) form a cover of M, 

(3) the multiplicity of the cover by balls B(£k,p) is not greater than 

Any set of points M p = which has properties (l)-(3) will be called a 
metric p-lattice. 

Let L be an elliptic second order differential operator on M, which is 
self-adjoint and positive semi-definite in the space L2(M) constructed with 
respect to Riemannian measure. Such operator has a discrete spectrum 
< Ai < A2 < .... which goes to infinity and does not have accumulation 
points. Let {uj} be an orthonormal system of eigenvectors of L, which is 
complete in L2(M). 

For a given uj > the notation E W (L) will be used for the span of all 
eigenvectors Uj that correspond to eigenvalues not greater than uj. 

Now we are going to prove existence of cubature formulas which are exact 
on E W (L), and have positive coefficients of the "right" size. 

The following exact cubature formula was established in [TT], [25j . 

Theorem 15. There exists a positive constant C, such that if 

P = C{uj + 1)- 1 I\ 

then for any p-lattice M p = there exist strictly positive coefficients 

p,£ k > 0, 6 M p , for which the following equality holds for all functions 
in E W (L): 

(23) / fdx= £ /^/fe)- 

Moreover, there exists constants c±, ci, such that the following inequalities 
hold: 

(24) cip n < < c 2 p n , n = dim M. 

9.2. On the product of eigenfunctions of the Casimir operator C 
on compact homogeneous manifolds. A homogeneous compact mani- 
fold M is a C°°-compact manifold on which a compact Lie group G acts 
transitively In this case M is necessarily of the form G/K, where K is a 
closed subgroup of G. The notation L2(M), is used for the usual Hilbert 
spaces L2(M, cZ£), where d£ is an invariant measure. 

Every element X of the (real) Lie algebra of G generates a vector field 
on M, which we will denote by the same letter X. Namely, for a smooth 
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function / on M one has 

t->0 t 

for every £ G M. In the future we will consider on M only such vector 
fields. The translations along integral curves of such vector fields X on 
M can be identified with a one-parameter group of diffeomorphisms of M, 
which is usually denoted as exptX, — oo < t < oo. At the same time, the 
one-parameter group exptX, — oo < t < oo, can be treated as a strongly 
continuous one-parameter group of operators acting on the space L2(M). 
These operators act on functions according to the formula 

/ -> f(exptX ■ £), t G R, / € L 2 (M), x G M. 

The generator of this one-parameter group will be denoted by Dx, and the 
group itself will be denoted by 

e tDx f(0 = f(exptX ■ £), t G R, / G L 2 (M), £ G M. 

According to the general theory of one-parameter groups in Banach spaces, 
the operator Dx is a closed operator on every L 2 (M). 

If g is the Lie algebra of a compact Lie group G then ([12], Ch. II,) it is 
a direct sum g = a + [g, g], where a is the center of g, and [g, g] is a semi- 
simple algebra. Let Q be a positive-definite quadratic form on g which, on 
[g, g], is opposite to the Killing form. Let X\, ...,Xd be a basis of g, which 
is orthonormal with respect to Q. Since the form Q is Ad(G)-invariant, the 
operator 

-X\ -X\- ... -Xj, d = dim G, 

is a bi-invariant operator on G. This implies in particular that the corre- 
sponding operator on L 2 (M) 

(25) C = -D\ -D\- ... - Dj, Dj = D Xj , d = dim G, 

commutes with all operators Dj = Dxj ■ 

This elliptic second order differential operator C is usually called the 
Laplace operator. In the case of a compact semi-simple Lie group, or a 
compact symmetric space of rank one the operator C is, or is proportional 
to, the Laplace-Beltrami operator of an invariant metric on M. 

The following theorem was proved in [TT], [25]. 

Theorem 16. If M = G/K is a compact homogeneous manifold and C is 
defined as in 125\) . then for any f and g belonging to E w (£), their product 
fg belongs to 'E^duj(C), where d is the dimension of the group G. 
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